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RESPONSE  OF  AN.  EIASTIC  DISK  TO  IMPACT  AND  MOVING  LOADS 


by 

A-  Cemal  Eringen  * 

Purdue  Unive  ra ity 

ABSTRACT 

With  the  use  of  FourLer  transforms  a class  of  elasto- 
dynamic  problems  concerning  disks  have  been  solved.  The  disk 
is  subjected  to  various  types  of  dynamic  loadings  at  the  rim. 
The  case  of  impact  and  the  moving  loads  are  studied  in  detail. 


(*}  Associate  Professor  of  Division  of  Engineering  Sciences. 


1.  Introduction 


Cylindrical  roller  bearings  in  high  speed  mechanism  are  subject 
to  dynamic  loadings.  Yet  the  usual  design  procedure  is  based  on  Hertz 
formulae  which  are  the  result  of  elasto-static  considerations.  In  many 
other  instances,  gears,  rollers,  or  disks  are  subject  to  impact  or  mov- 
ing loads.  If  we  neglect  the  coriolis  terms  we  can  also  bring  the  roll- 
ing disks  on  contact  into  the  category  of  disks,  subject  to  moving  loads. 
Thus  the  aim  of  the  present  paper  is  to  obtain  the  solution  to  this 
class  of  elasto-dynamic  problems  concerning  the  disk.  The  dynamic  load 
Is  applied  to  the  rim  of  the  disk.  Two  normal  concentrated  dynamic 
loads  at  the  two  ends  of  a diameter  moving  or  otherwise  are  special 
cases. 

Some  solutions  of  free  oscillation  of  cylinders  are  known  since 
, the  time  of  Pochhammer  /~1_7»  and  later  Pickett  /~2_7,  J.  Uindlin  /~3_7> 
T.  Ghosh  Similarly,  the  problem  of  rotating  disks  has  attracted 

attention  of  many  authors  (see,  for  instance,  Lamb  and  Southwell  /~5_7> 
Timoshenko  and  Goodier  /~6_7»  Love  /”7_7  )•  It.  seems,  however,  that  the 
forced  oscillation  problems  concerning  disks  and  cylinders  have  escaped 
the  attention  of  authors,  excepting  a paper  by  J.  Uindlin  /~8_7,  which 
consists  of  generalities  in  a related  problem  to  ours. 

The  present  method  is  applicable  to  ring  problems  and  to  the 
plane  elasto-dynamic  problems  concerning  circular  holes.  These  prob- 
lems will  be  treated  in  later  papers. 


2.  Formulation  of  the  Problem 


Equations  of  motion  of  plane  homogeneous  isotropic  media  in  terms 
of  plane  polar  coordinates  r and  0 and  the  time  t are  / 7_7:  (See  Fig.  1) 
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Fig.  1.  Circular  Disk 
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where  u (r,  0,  t)  and  v (r,  0,  t)  are  components  of  the  displacement 
vector,  X and  are  the  Lame^  constants,  and  Y is  the  mass  density 
per  unit  volume.  Dilatation  A and  rotation  cO  are  related  to  u,  v 
by: 


r A ■ (r  u) 


v,0  , 2 r tJ  - ( r v)^  - u^  (2) 


Subscripts  after  a conma  represent  differentiation,  i.e.,  u 
^ ^ 0 etc . .... 


Elimination  of  u and  v among  (1)  and  (2^  leads  to: 


\r\r  • ‘,86  ' ’ V /(  X • 2 jX ) 


OL^  r^cO  ..  - r (r^L>  ) * eJ 

2 »tt  ,r  ,r  ,0Q 
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> 0^2 


y //X 


(3) 


These  are  the  equations  of  dilatatlonal  and  rotational  waves. 

Components  0"  , ff*  , 0^^  of  the  stress  tensor  are  given  by: 
rr  rO 

<r„  ’ X a*  2^1  - |X  r-1  ^ r(v/r)^^ 

- X A . 2^r-l  (v^^.  u) 


(4) 


The  problem  is  to  solve  (1),  (2)  under  a given  and 

at  the  riffl  r ■ a of  the  disk: 

<r^  (a,  0,  t)  - fi;  («,  t)  , (Tg  (a,  0,  t)  - r (©,  t)  (5) 

subject  to  the  condition  that  these  surface  tractions  are  in  equilibrium 
at  each  Instant. 


3.  The  Solution 

The  periodic  solution  of  (3)  with  respect  to  6 is  obtained 

to  be 

_ oo 

4 ■ ^ (A,  sin  n 0 * cos  n 6)  Z_  ( P, ) , 

^ In  2n  ' n ' / 1''  ’ 


n«o 


^ cos  n 0 - sin  n ©)  ( p^)  , 


(6) 


n“o 
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where  Zq(^)  ie  the  cylinder  function  tS’J*  ^jn'  ®jn'  ^Jn  ^Jn 


are 
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constants  of  integration;  and  the  barred  quantities  represent  the 
Fourier  transforms,  i.e. 


F it) 


f e^^^  F (t)  dt  , i « l/^ 

•'-oo 


(7) 


The  inversion  formula  for  (7)  is 

lOO 


^ ’ i;;  f ^ (r)  dr 

2tt  J_oo 


(8) 


Substituting  (6)  into  (2)  after  taking  the  Fourier  transforms 
of  (2)  and  solving  the  resulting  partial  differential  equations  we 
obtain: 

“ ■ H “in''-)  sin  n VJ  ♦ ^2n^’^^  ^ 


oo 
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n»o 

oo 


Vin(r)  cos  n « - ^ 


n-o 


2n' 


where 


-r 


-1 


V--’  ■ *jn  •'f‘l 


-1 


where  prime  represents  differentiation. 

Combining  (4)»  (6),  and  (9)  we  obtain 
oo 

• ^^*ln  V''*’’  * “in  " ® '^*2n"ln< 

oo 

“ln<«''  * “in  " ® '^‘anV*’ 

(cont'd. ) 


(9) 


(10) 


(11) 
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(11,  cont'd.) 
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where: 

«ln  (■•Z') 

(rr) 
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4,  Dynamic  tractions  applied  to  the  rim  of  a dlak 

In  the  case  of  a disk,  the  stress  and  deformation  components 

must  be  finite  at  r ■ 0.  Hence  ■ 0.  Without  loss  of  generality 

we  also  take  - 1.  This  means  in  all  of  our  formulas  we  must  re- 

Jn 

place  by  J^. 

We  use  boundary  conditions  (5)  to  determine  the  constants  Aj^ 

and  B.  . The  Fourier's  theorem  thus  leads  to: 
jn 

*ln  “n  Ci2n  - »2„ 
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where: 
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<r  (0,r)  cos  n 0 de 
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Various  special  cases  are  of  Interest: 


(a)  Zero  surface  shear;  Z*  ■ Z"  • 0 

" OC  08 

(b)  Normal  traction  with  central  symmetry  and  (a) ; - 0 

(c)  (b)  with  constant  amplitude  over  o - Q ^ oC  : 


(15)  and  (T^  (©,r)  ’ it)  J 


0 ^ 
n - 0^  - © ^ TT 


Hence 
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(d)  Impact  load  and  (cis 
(16)  and 


lim  a «■  (r)  a OC  - p (r) 

0 0 
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Hence 
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“■oc  ■ f. 

(e)  Impulalve  concentrated  load  and  (d)t 

(18)  and  (t)  • **o  where  5”  (t)  is  the  Dirac  delta 


function  defined  by: 

(t)  « ^ oo 
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In  this  case« 
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where  P^  is  the  amplitude  of  the  concentrated  load. 


(19) 


5.  Moving  load 

(a)  Moving  normal  and  tangential  stresses: 

Moviiig  loads  can  be  represented  by 

<r  (0,  tl  • (r(g-^t)  , r (o,  t)  - r(g_^t) 

o A o 0 ox 

Fourier  Transforms  of  these  are: 


(20) 


^ (0,  r ) - 2 a^(«'A»-)  , r^(0,r)-  2e^^®/-^l  8^(r/A^) 


(21) 


s (r//l)  - <5l(0)e  0 (in  general) 

o ‘■‘^-00  ° 

~ <r(0)co8  (T0/a)d0  when  d^(-0)  - 6^(0) 

where  follows  from  8^(  r'/zi)  by  replacing  U7.  and  <r(0) 

by  and  ^(0)  respectively.  Haice: 


(22) 


(23) 


- (2r/n^\)  Z"(r/^)^  - s^crAz-) 

? - (2  n/n^)  r^XjU^f  - (e^^^^/-^-l)  s (r/ri) 

ua  o 

Quantities  T and  ^ are  obtained  from  (23)  by  writing-/?  and 
oc  os  X 

{t  / in  place  of and  ®q  (7*  A^- ) respectively. 

(b)  Moving  periodic  loads; 

oo 

<r  (e  t)  - (Q  /naj  p cos  n(0  -./2t;  + q sin  n(Q--/J.t) 

° n“o 

(24) 
oo 

^ (0  - -/2^t)  - (Qo/Tia)^  r^  cos  n(0  --/2^t)+s  sin  n(0-  -^^j^t) 

n“o  ^ 

After  taking  Fourier  transforms  of  <7"  and  "Z"  > w®  substitute  into 

o o 

(13)  to  obtain  <r  , , f'  . This  gives; 

OC  OS 

® ^Pn  ■ Z (/In  - T ) ♦(  p ♦ iq  ) S’  (-J7n  -D 

a 6"  /Q  - (ifp  ♦ q ) S’  (iT.n  -T)  * (-ip  ♦ q ) S’  (-.nn  -t ) 
o®  o n n n n 

Quantities  t'  and  X follow  from  (25)  by  witing  r„  and 
oc  os  n 

place  of  p^  and  q^  . In  obtaining  (25)  we  used  the  formal  relation* 

2n  S’  (u)  - f dt 

•'-oo 

Components  of  displacements  and  stress  tensor  csui  now  be  obtained  by 
combining  (9),  (11) > (13)  and  (25)  and  taking  inverse  Fourier  trans- 
forms . Thus 


(25) 

s in 
n 


♦ This  is  Justified  in  the  sense  of  a distribution  function  /~11  7. 


-(2rT  a W/Q^r)  u r^n)/”p  cos  n (0-^t)+q^sin  n(0-^t)_7 

° n*o 


/ O \ 

\rjT.  n)^r^sin  n (©-i2^t)  + 8^cos  n(0-j7^t)_7 


- (2n  a l/^  r)  v v^^\r^n)/^p^sin  n(Q--/lt)+q^cos  n(9-^t)_7 

^ ° n=o 

)/?r  ^cos  n(Q-/1^t)-SjjSin  n(0-/2^t)_7 


(n  a/Q  )0'  - X 6^^\rAn)/"p„cos  n(0-At)+q  sin  n(«-/it)_7 

' o rx*  • # xn  n n 

n=o 

+ d^^^(rn,n)^r  sin  n(e-y2  t)+s  cos  n(9--«7it)_7 
In  i n 1 “ 

OO  (1  ) _ 

(n  a/Q  )fll^  - ^ T (ri:2n)^p  sin  n(9-^t)+q  cos  n(9-At)_7 
' o r«  n n “ 


+ {TJi.n)/Pr  ^cos  n(0-j(^t)  - s^sin  n(9-^t)_7 

(n  a/QQ)<5^  6'2^\r^n)/’p^cos  n(0-j7t)+  q^sin  n(9-ilt)_7 

•*■  6^^  (r^n)/-r  sin  n(0-yi,t)  + s cos  n(9-/i  t)_7 
‘ 2n  n T.  n 1 


where 


■ (a^n)J7  ruin)  (oc'^r^n) 

- 2n(<X2r^n)  (o/^r^n)  (axin)_/ 

-1  _2 

" Z~Dfj  (a^n)_7  Z~n  (o^^r^n)  ((Xj^ri2n)  S2jjaxi.n) 

- 2 (cv^r^n)  ^ ( o^2rxin)  (a^n)_7  (27) 

-1 

- /Djj  (a^n)_7  ^j^^(r^n)  S^^Ca^n)  -N^j^Cr^n)  S^^(a^n)_7 
-1 

(ax2n)_7  /S^^(rj^n)  S^^(a^n)  -S^^(rj^n)  S^(a^n)_7 
■(j^\r^n)  - (aAn)_7  S^^Ca/in)  -T^j^Cr^rzn)  S^^(axin)_7 

where 

Djj(a^n)  - Nj^^(auzn)  S2jj(a^n)  - N^^Ca^n)  S^(a^n)  (28) 

Functions  u^^^,  \ arid  are  obtained  from  the 

n n In  n 2n 

corresponding  ones  with  superscripts  (l)  above  by  replacing^  by  ^ ^ 
and  S2jj(a^n)  and  S^jj(a^n)  by  N2jj(a-n^n)  and  respec- 

tively, except  in  D^(aAn),  where  we  replace  hy . 


u^^\r^n) 

n 


v|j^\riin) 


ff  (rjan) 


/ 


(c)  Two  diametrically  oppoalte.  moving,  concentrated  loads 
(Fig.  2). 

This  is  a case  of 
technical  Importance. 

In  this  case  we  have 


7'  • 0 hence 
o 

“ a„  - 0.  The  con- 
n n 

centra ted  loads  can 
formally  be  represent- 
ed by 

6^(9- yzt)  - (Q^n  a)  /"/(0-JT.t)  + (T  (n  - 9 * -n.t)J 

oo 

z 

n-0 


(29) 


(Q^n  a)  JZ  Pn  \ -^t) 


where  is  the  amplitude  of  each  of  the  concentrated  radial  load. 

From  (29)  Fourier  coefficients  p^  and  ^ are  calculated  to  be: 
2/n  for  n even 


for  n odd 


^n  ' 0 


(30) 


Hence  (26)  together  with  (29)  and  (30)  and  r = s - 0 gives 

n n 

the  displacement  and  stress  components.  Below  we  give  displacement  com- 
ponents. The  rest  is  obtained  in  an  obvious  manner. 

.(1)/ 


- (n^  a li/Q  r)  u - ^ u ^ (ry2n)  cos  n (9  - JXt) 

0,2,4,...  ” 

2 (1) 

(n  a ^/Qq  r)  V - y , v^  (r^n)  sin  .i  (9  ~yit) 

0,2,4, ... 


(31) 


where  ouid  v^^^ 

n n 


are  given  by  (26). 


Concentrated  moving  shear  loads  and  other  types  of  load  combina- 
tions nay  easily  be  obtained  from  (24)  and  (26). 

6.  Computation  and  Discussion 

Computations  have  been  carried  out  to  determine  the  roots  of 
D^(aj2n)  ■ 0 given  by  (28).  The  roots  of  this  equation  give  the 
resonance  speed  for  the  concentrated  load. 

The  ratio  of  the  resonance  speed  C^  • a jOl  of  the  concentrated 
load  to  the  dilatational  wave  velocity  C^^  • ^1~^  been  solved  from 
- 0 for  a steel  cylinder  with  E * 30  x 10^  psi  and  V * 0.3*  First 

four  roots  of  ■ 0 for  n ” 1,  2,  3>  4 are  listed  on  the  table  given 

* 

below  . 

The  computation  is  carried  out  on  l.B.U.  Card  Program  Calculator 
machine  and  are  correct  up  to  two  decimal  places. 

The  values  of  n C^C-^  ~ n (obtained  by  calculating  the  roots 

of  Djj  (a-an)  - 0).  Stool;  E = 30  x 10^  , V - 0.3  . 


n 

= 0 

D2  = 0 

Dj-O 

°4  = 

1 

1.51 

1.26 

1.94 

2.52 

2 

3.43 

2.35 

3.21 

4.08 

3 

3.80 

4.22 

4.95 

5.67 

4 

5.33 

5. 06 

6.22 

7.19 

* For  n “ 0 we  take  the  limit  of  the  functions  in  (27)  and  find  finite 


values 


From  this  table  it  is  seen  that  the  smallest  is  obtained 

for  n-2.  This  «lue  is  CjyC^-  1.26/2  - .63.  We  therefore  expect  a 
critical  speed  for  the  moving  load  in  the  neighborhood  of  O.63  which 
will  create  resonance  in  the  disk. 

The  Rayleigh  surface  w^ve  velocity  for  1)  “ .29  is  0.9258  times 
the  velocity  of  shear  waves  or  .503  times  the  velocity  of  dilational 
waves.  Thus  0.63  represents  a velocity  between  the  Rayleigh  surface 
wave  velocity  and  the  shear  wave  velocity.  Further  we  notice  that  this 
velocity  is  minirmm  for  n-2  rather  than  n-1.  Examining  (31)  we  found 
that  we  have  only  even  terms.  Hence  diametrically  opposite  moving  loads 
give  a smaller  critical  speed  than  the  case  of  a single  load.  In  the 
former  case  the  fundamental  mode  n«0  represents  a tiniform  lateral  exten- 
sion which  is  not  dependent  on  time.  Thus  it  is  static  in  nature  and 
has  no  resonance  frequenqy  associated  with  it.  The  second  mode  is  of 
cos  2 0 — type  and  gives  the  minimum  critical  speed  mentioned  above. 

It  is  interesting  to  know  that  the  minimum  critical  speed  in  the  case 
of  a flat  semi-infinite  plate  (the  Rayleigh  surface  wave  velocity)  is 
less  than  that  of  the  curved  surface  ( the  disk  in  the  present  case). 
With  this  point  in  mind  pertiaps  one  can  classify  the  curved  surfaces 
from  the  elasto-dynamic  point  of  view.  This  process  will,  no  doubt, 
present  mathematical  difficixlties. 

The  present  results,  of  course,  must  be  accepted  tentatively 
until  experimental  verification. 
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